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ABSTRACT 

We present the results of a search for high energy neutrinos with the Baikal 
underwater Cherenkov detector NT-200. An upper limit on the (i^e + J^e) diffuse 
flux of E'^^^{E) < (1.3^ 1.9) • 10~^ cm~^ sr~^ GeV within a neutrino energy 
range 10^ 10^ GeV is obtained, assuming an E^"^ behaviour of the neutrino 
spectrum and flavor ratio (i^e + i^e) ■ {vp, + i7^)=l:2. 

1. Introduction 

The Baikal Neutrino Telescope is deployed in Lake Baikal, Siberia, 3.6 km from 
shore at a depth of 1.1 km. The optical properties of Lake Baikal deep water are 




Figure 1: Schematic view of the Baikal neutrino telescope. 



characterized by an absorption length of 20-=- 25 m, a scattering length of 20-^70 m 
and a strongly anisotropic scattering function f{6) with mean cosine of scattering 
angle cos(6')=0.85-r0.95. NT-200, the medium-term goal of the collaboration was 
put into operation on April 6th, 1998 and consists of 192 optical modules (OMs). An 
umbrella-like frame (see Fig.l) carries 8 strings, each with 24 pairwise arranged OMs. 
Three underwater electrical cables and one optical cable connect the detector with 
the shore station. 

The OMs are grouped in pairs along the strings. They contain 37-cm diame- 
ter QUASAR - photomultipliers (PMs) which have been developed specially for our 
project I-*. The two PMs of a pair are switched in coincidence in order to suppress 
background from bioluminescence and PM noise. A pair defines a channel. 

A muon-trigger is formed by the requirement of >4 hits (with hit referring to a 
channel) within 500 ns. For such events, amplitude and time of all fired channels are 
digitized and sent to shore. Full space-angular track reconstruction is possible for 
muon induced events with >6 hits at >3 strings {off-line trigger 6/3). A separate 
monopole trigger system searches for clusters of sequential hits in individual channels 
which are characteristic for the passage of slowly moving, bright objects like GUT 
monopoles. 

Table 1 illustrates A'^T-^Of operation efficiency since April 1998 till February 2001. 
Shown are the data taking time T, the fraction T/Ttot of total detector operating time 
Ttot when data have been taken, the number of accumulated events Ng^, the number of 
events allowing full space-angular track reconstruction Ne^(6/3), the average fraction 
of working channels l^op/'^tot as well as the data taking time T(>85%) with >85% 
channels in operation. 



Table 1: NT- 200 - operation efficiency 
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T 

(days) 


T/Ttot 


(10^ events) 


Ne.(6/3) 

(10^ events) 




T(>85%) 
(days) 
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73% 


167 


57 


71% 
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99-00 
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75% 


191 


65 


76% 


61 


00-01 
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79% 


233 


82 


81% 
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Using the data accumulated during first 234 livetime days between April 1998 
and February 1999, 35 neutrino induced upward muons have been reconstructed i*. 
Although in a good agreement with MC expectation this number is by a factor 3 
lower than predicted for the fully operational NT- 200. The reason is that, due to 
unstable operation of electronics, in average only 50 - 70 channels have taken data 
during 1998. This is in contrast to 1999 and 2000 data taking where the stability 
had been improved. Ten events within a 30 degree half angle cone around nadir have 
been selected to set limits on the excess of the muon flux due to WIMP annihilation 
in the center of the Earth. Also a new limit on the flux of fast monopoles has been 
obtained & . 

Here we present preliminary results from a search for high energy neutrinos 
(E;, >10 TeV) with NT-200 obtained from the analysis of the entire 1998 data set. 

2. A search for high energy neutrinos 

The main goal of large underwater neutrino telescopes is the search for extrater- 
restrial high energy neutrinos. Detection volume and detection area of such detectors 
depend on the transparency of the medium (water or ice) and the power of the source 
of Cherenkov radiation (high energy shower or muon), and may significantly exceed 
the geometrical one. 

2.1. Search strategy 

The used search strategy for high energy neutrinos relies on the detection of the 
Cherenkov light emitted by the electro-magnetic and (or) hadronic particle cascades 
and high energy muons produced at the neutrino interaction vertex in a large vol- 
ume around the neutrino telescope. Earlier, a similar strategy has been used by 
DUMANdS), AMANDA i) and BAIKAL & collaborations to obtain upper limits 
on the diffuse flux of high energy neutrinos with relatively small detectors (SPS, 
AMANDA-A and NT-96, respectively). Although the limits implied by these obser- 
vations are at least one order of magnitude higher then the model independent upper 
limit derived from the energy density of the diffuse X- and gamma-radiation &, these 
results illustrate the power of used search strategy. 



Neutrinos produce showers through CC-interactions 



+N ^ /-(/+) + hadrons, (1) 

through NC-interactions 

z/,(z7,) + z/,(z7i) + hadrons, (2) 

where I = e, fx or t, and through resonance production iSB) 

z7g + e~ — > W~ — > anything, (3) 

with the resonant neutrino energy Eq = M^/2me = 6.3 ■ lO^GeV and cross section 
5.02 ■ lO-^^cm^. 

We select events with high multiphcity of hit channels N^it corresponding to bright 
cascades. The volume considered for generation of cascades is essentially below the 
geometrical volume of NT-200. A cut is applied which accepts only time patterns 
corresponding to upward traveling light signals (see below). This cut rejects most 
events from brems-cascades produced by downward going muons since the majority of 
muons is close to the vertical; they would cross the detector and generate a downward 
time pattern. Only the fewer muons with large zenith angles may escape detection and 
illuminate the array by their proper Cherenkov radiation or via bright cascades from 
below the detector. These events then have to be rejected by a stringent multiplicity 
cut. 

The used strategy is very effective for z/g and detection since the main fraction 
of neutrino energy would be transfered to electro-magnetic or /and hadronic cascades 
due to CC-interactions. It is less effective for z/^ detection since the main part of 
neutrino energy would be escaped by energetic muon from detection volume. For 
search preliminary result has been presented by the AMANDA experiment % . 



2.2. Data 

Within the 234 days of the detector livetime, 1.67 • 10^ events with N^n > 4 have 
been selected. For this analysis we used events with N^j^ >10. The time difference 
between any two hit channels on the same string was required to obey the condition: 



(t,-t,) > -10 ns, («<j), (4) 

where tj, tj are the arrival times at channels i, j and the numbering of channels rises 
from top to bottom along the string. 

Since April 1998 till February 1999, NT-200 operated in 3 main configurations 
with 77, 60 and 49 working channels, respectively. Data taking time T, the number 
of selected events Net, which survive cut (4) and the maximum hit multiplicity ^^^^ 
of these events, are shown in Table 2. 
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Figure 2: Distribution of hit channels multiplicity; dots - experiment, hatched boxes - expectation 
from brems and hadronic showers produced by atmospheric muons. 



Table 2: NT-200 configurations in 1998. 



Configuration 


Nop 


T 

(days) 


Ne^ 




Nthr 


1 


77 


57.9 


63540 


45 


50 


2 


60 


145.8 


83319 


37 


39 


3 


49 


30.9 


8719 


24 


26 



Fig. 2 shows the hit multiphcity distribution of selected (dots) as well as the ex- 
pected one from background high energy brems and hadronic showers produced by 
atmospheric muons (hatched boxes). The experimental distribution is consistent with 
the theoretical expectation for ^Mt >18. For lower N^jj values the contribution of 
atmospheric muons close to horizon as well as low energy showers from e^e~ pair pro- 
duction become important. The highest multiplicity of hit channels experimentally 
observed is NJ^^^ = 45 (one event). No statistically significant excess over expectation 
from atmospheric muon induced showers has been observed for each of the 3 detec- 
tor configurations. The detection efficiency of NT-200 for events with Nhu > Nf^^^ 
had been analysed by applying several less stringent cuts. It was shown that the 
experimental distributions are consistent with expected ones from atmospheric 
muons. 

Since no events with Nhu > NfJ}^^ are found in our data we can derive upper 
limits on the flux of high energy neutrinos which would produce events with 
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Figure 3: The normalized hit muhiphcity distributions of events which would be produced by the 
fluxes and survive the selection criterion (4). Solid, dashed and dotted curves correspond to 7=1.5, 
2, 2.5, respectively. Also shown is the normalized ^hit distribution of events from atmospheric muon 
induced showers (thick line). 



where the chosen values of ^thr for the 3 detector configurations are given in Table 2. 

2.3. MC- simulations 

Given an isotropic diffuse high energy neutrino flux with power law energy spec- 
trum with spectral index 7, the number of expected events during observation time 
T reads 

N, = ^JdnJ dEVeffin, E)J2j dE,E;mAPH,o^ exp{-liQ)/hot) (6) 

k 

where Ei, is the neutrino energy, E - the energy transferred to a shower, A,y - normal- 
ization coefficient of neutrino fiux and Veff{^l,E) - detection volume 0. The index 
u indicates neutrino types (z/ = u^, i7^, Ue, Ue) and k indicates the summation over 
CC and NC interactions. A^^ is the Avogadro number. Cross sections 0© da^k/dE 
correspond to processes (1) and (2). The neutrino absorption in the Earth has been 
taken into account with a suppression factor exp(— Z(f2)/Ztot), where 1{VL) is the neu- 
trino path length through the Earth in direction VL and /^i = PEarth [o'cc + o"Arc) 
according to 0©. 

In Fig. 3 we show the normalized predictions of '^hu distributions of events which 
survive cut (4) and would be induced by electron neutrino fiuxes with 7=1.5, 2, 2.5. 
The normalized distribution of events from atmospheric muon induced showers, 
which has strongly steeper behaviour, is also presented. 

^In a case of Vfj_ CC- interaction in water the detector response to hadronic shower as well as to 
the high energy muon has been taken into account. 
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Figure 4: The fractions of expected events induced by diffuse Ve fluxes witli spectral indexes 7=1, 
1.5, 2, 2.5 within energy range 10 TeV<E^ <E. 
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Figure 5: Detection volume of Ve produced events which survive cuts (4)-(5) (upper curve). The 
detection volume folded with the neutrino absorption probability in the Earth (lower curve) is also 
shown. 

The neutrino detection energy range of NT-200 which contains, for instance, 90% 
of expected events, depends on the value of spectral index 7. Normalized energy 
distributions of expected events induced by electron neutrinos within energy range 
10 TeV E (E<10^ TeV) and different 7 are presented in Fig.4. Assuming 7 ~2 as 
typically expected for Fermi acceleration, 90% of expected events would be produced 
by neutrinos from the energy range 20 ^ 10^ TeV with the mean energy around 200 
TeV. 

The detection volume for neutrino produced events which fulfill conditions (4)- 
(5) was calculated as a function of neutrino energy and zenith angle 9. The energy 
dependence of the detection volume for isotropic flux with 7 —2 is shown in Fig. 5. 
Also shown is the detection volume folded with the neutrino absorption probability 
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Figure 6: Coordinates of interaction vertex for events which fulfiU conditions (4) (dots) and (4)- (5) 
(rectangles). 

in the Earth. The value of Veff rises from 2-10^ for 10 TeV up to 6-10^ for 10^ 
TeV and significantly exceeds the geometrical volume ~ 10^ m'^ of NT- 200. This 
is due to the low light scattering and the preserved light fronts from Cherenkov waves 
originating far outside the geometrical volume. Although the detection volume has 
been calculated without taking into account light scattering in the water, estimations 
show that a scattering with 1^=20 m and cos(6')=0.88 (conservative values for Lake 
Baikal water) would cause <30% decrease of Veff for <10^TeV. 

Fig. 6 illustrates the difference between Veff and Vg. Shown here are the coor- 
dinates of neutrino interaction vertex for events which survive cuts (4) (dots) and 
(4)- (5) (rectangles). 

2.4- The limits on the high energy neutrino fluxes 

The shape of the neutrino spectrum was assumed to behave like E""^ and flavor 
ratio (z/e + z/g) : {^fi + ^v) = 1:2 due to photo-meson production of vr"*" followed by 
the decay Tr~^ ^ fj,~^ + ^ + + i^fj, + i'^ for extraterrestrial sources. 

Comparing the expected number of events fulfilling (4)- (5) with the upper limit 
on the actual number of events, 2.4 for 90% CL. © we obtain the upper limit on 
the diffuse (z/g + z7g) flux. The combined upper limit obtained with Baikal neutrino 
telescopes NT-200 (234 days) and NT-96 i (70 days) is: 

^iu.+u,)E^ < (1.3 - 1.9) • lO-^cm-^s-isr-^GeV, (7) 

where the upper value refers to the conservative limit on light scattering in the Baikal 
water. 
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Figure 7: The upper limits on the diffuse neutrino fluxes, obtained by BAIKAL (this work: for 
flux as well as for a model independent flux at resonant energy 6.3-lO^GeV (rectangle)), 
AMANDA and FREJUS (triangle) experiments as well as the atmospheric conventional neutrino 
fluxes from horizontal and vertical directions (upper and lower curves, respectively) and atmospheric 
prompt neutrino flux (curve labeled Vpr)- The limits and predictions from different models of high 
energy neutrino sources are also shown. Curves labeled 'B', 'MPR' and 'WB' represent the model 
independent limit derived by Berezinsky and the upper bounds derived by Mannheim et al. as well 
as by Waxman and Bahcall. Curves labeled as 'SS' and 'P' show predicted fluxes from quasar cores 
and Protheroe respectively. Curves labeled 'M(GRB)' and 'WB(GRB)' represent predicted fluxes 
from GRBs derived by Mannheim and Waxman and Bahcall. Curve labeled 'TD' - prediction for 
neutrino flux from topological defects due to specific top-down scenario BHSl. 



Fig. 7 shows the upper hmits on the isotropic diffuse high energy neutrino ffuxes 
obtained by BAIKAL (this work), AMANDA i and FREJUS 0) (triangle) experi- 
ments as well as the atmospheric conventional neutrino ffuxes from horizontal and 
vertical directions (upper and lower curves, respectively) and atmospheric prompt 
neutrino ffux 0^ (curve labeled Vp,). Also shown is the model independent upper 
limit on the diffuse high energy neutrino ffux obtained by V. Berezinsky % (curve la- 
beled 'B') with the energy density of the diffuse X- and gamma-radiation cu^ < 2 • 10~^ 
eV cm~^ (as follows from EGRET data©), and predictions for diffuse neutrino fluxes 
from Stecker and Salamon model 0* (the sum of the quasar core and blazar jet con- 
tribution, curve labeled 'SS') and Protheroe model El* (from equal contribution from 
pp and p7 interactions in AGN jets, curve labeled 'P'). Curves labeled 'MPR' and 
'WB' show the upper bounds obtained by Mannheim, Protheroe and Rachen for 
optically thick {jn^ > 1) and optically thin (r„,^ < 1) sources as well as the upper 
bound obtained by Waxman and Bahcall 0^ for optically thin sources respectively. 
Curves labeled 'M(GRB)' and 'WB(GRB)' present the upper bounds for diffuse neu- 
trino flux from GRBs derived by Mannheim 0^ and Waxman and Bahcall il^. Curve 
labeled 'TD' shows prediction for neutrino flux from topological defects due to specific 
top-down scenario BHSl 0). 

For the resonant process (3) the event number is given by: 



(M,„+2r,„)V2me 

AT,-, =TjdnJ dEVejfiQ, E) J dE,^,-XE,)^NApH,o^ (8) 

= 80.22GeV, = 2.08GeV. 

Our combined 90% C.L. limit obtained with NT-200 and NT-96 i at the W - 
resonance energy is: 

^ < (1.4 - 1.9) X lO-^'cm-'s-^sr-^GeV^^ (9) 

Cl-C/p 

and is also depicted in Fig. 7 (rectangle). 

3. Conclusion 

The deep underwater neutrino telescope NT-200 in Lake Baikal is taking data 
since April 1998. Due to the high water transparency and low light scattering, the 
detection volume of NT-200 for high energy z/g and z/^ detection is several megatons 
and exceeds the geometrical volume by factor of about 50 for highest energies. This 
permits a search for diffuse neutrino fluxes from extraterrestrial sources on the level of 
theoretical predictions. The upper limits (7), (9) obtained for the diffuse (z/g + t'e) 
flux and the model independent Ve flux at resonant energy 6.3-lO^GeV are the most 
stringent ones at present. We expect that the analysis of 3 years data taken with 
NT-200 would allow us to reach a sensitivity of ^^E"^ 6 • 10^''cm~^s~^sr^^GeV. 
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